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1. Introduction 

Real exponential functionals of Brownian motion have been the subject of much activity 
in probability theory [1] . They have found many applications to problems ranging from 
finance to physics [2]. To be more specific, if B{t) is Brownian motion (a Wiener process 
with S(0) =0 and {B{tY) — t), the following stochastic variable 

Xit) = /*e-*+^^Wds, (1.1) 

^0 

where is a real coupling constant, shows up in one guise or another in various models 
of disordered systems [3, 4]. The above random integral can be shown to represent the 
solution of the following Langevin equation with multiplicative noise: 

^^l-X{t)+gXmt\ (1.2) 

with Tjit) = dB(t)/dt being a zero-mean Gaussian white noise, normalized as 

{v{sMt)) = 5{t~s). (1.3) 

The derivation of the solution to (1.2) in the form (1.1), including the required stochastic 
calculus prescription (Stratonovich), will be reviewed in Section 2. As an illustration of 
the versatility of the situations in which the random variable X{t) appears, the Langevin 
equation (1.2) may serve to study the effect of multiplicative noise on a deterministic 
fixed point. Indeed, the linear force 1 — X{t) stabilizes X{t) at the fixed point X — 1 
in the absence of noise. The solution X{t) however keeps on fluctuating forever because 
of the multiplicative noise term gX{t)r]{t). In the long-time limit, what emerges out of 
this interplay between the deterministic restoring force and the fluctuating noise term 
is a non-trivial distribution for the random variable 

X = / e-*+^^(*)di, (1.4) 
Jo 

representing the stationary solution of (1.2). The distribution of X = limt^oo-^(^) can 
be computed in a number of ways, and we shall review in Section 2 two useful methods 
to do so. Using either more elaborate path-integral methods or probabilistic identities, 
the full-time dependent distribution of X{t) can also be obtained, and features, among 
other things, an interesting continuous spectrum of relaxation rates [2] . 

A seemingly innocuous and somewhat natural generalization of (1.4) consists in 
analytically continuing the coupling constant g ^ \g. We thus obtain the following 
random integral, defining an imaginary exponential functional of Brownian motion: 

Z= \ e-*+^^^Wdt, (1.5) 

which now lives in the complex plane (actually in the unit disk). 

The main goal of the present article is to investigate the distribution of the 
functional Z. Throughout the following, instead of the coupling constant g^ we rather 
make use of the dimensionless diffusion constant 
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As hinted at above, even though a large number of works have been devoted to real 
exponential functionals of Brownian motion, much less is known about the distribution 
of complex functionals of Brownian motion such as (1.5). We found interesting from a 
conceptual viewpoint to tackle this problem, in particular in order to see if any of the 
methods which proved successful for the real case would extend to this complex- variable 
setting. Besides this, it turns out that complex stochastic processes have surfaced time 
and again in different scientific disciplines ranging from signal theory, where processes 
involving imaginary exponential functionals of Brownian motion occur in the study of 
phase noise [5], to quantum optics [6], in conjunction with the development of phase- 
space representations and of the associated formalism of quasi-probabilities, and finally 
to reaction-diffusion processes [7, 8] through the Doi-Peliti approach [9]. The latter 
topic has constituted our original thrust to embark on the study of (1.5). Let us now 
describe how the connection emerges. 

In the context of interacting particle systems such as reaction-diffusion systems, 
there exists a standard set of techniques, usually referred to as the Doi-Pehti 
formalism [9] (see [10] for a recent review, and the references therein), which allows 
to recast the master equation describing the evolution of these processes in terms of 
a field theory whose action involves a pair of conjugate fields. Without entering into 
much detail, provided certain technical conditions are met, the theory can in turn be 
transformed into a Langcvin equation for a single density field ^Pn{t), customarily dubbed 
so because its noisc-avcragc ((/9„(t)), coincides with the local mean particle number for 
the underlying reaction-diffusion process ((pn(t)))- Here the double brackets denote an 
average over the dynamics of the particles. In the case of the binary annihilation reaction 



where particles A diffuse by hopping with rate k on a hypercubic lattice in dimension d, 
and annihilate pairwise with rate A when they meet on a given lattice site n, one can 
show [7, 8, 9] that the stochastic density field (pn{t) obeys the following Langevin-Ito 
equation: 



with being the lattice Laplacian, such that, e.g., (V^</7)„ = (fin+i + fn-i — 2</?n 
in one dimension. It can be expected on physical grounds that the amplitude of the 

Gaussian noise term Cniv'n] vanishes when = 0. Indeed, the first two terms on the 
right-hand side of (1.8) respectively account for the diffusion of the particles and for the 
(mean-field) decay rate of the particle density due to pairwise annihilation. If the noise 
amplitude vanishes when = 0, there is no evolution at all in regions where = 0. 
The Doi-Peliti approach also yields the following expression for the correlator of the 
noise: 



where rin{t) is a normalized real Gaussian white noise. This explains why (1.8) 



A + A^(J}, 




dt 



/«(VV)n(i) -A(^„(i)' + Cn(i), 



(1.8) 
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is often referred to as an imaginary-noise equation. One would obtain the same 
Langevin equation using Gardiner's quasi-probability formalism, where the particle- 
number probability distribution is represented as a superposition of Poisson distributions 
with weights [8, 11]. The equivalence between the Doi-Peliti formalism and 
Gardiner's Poisson representation method has been demonstrated in general in [12]. 
There is no contradiction in either formalism, as soon as the auxiliary field (pnif) 
is complex- valued, provided one refrains from erroneously identifying (pn{t) with the 
(integer- valued) stochastic variable Pn{t)-i based on the sole equality between the mean 
values (</?„(i)) = {{pn{t))). In fact, the precise relationship between the distribution of 
(fin{t) and that of p„(t) is that the ordinary moments of (pnit) are equal to the factorial 
moments of Pn{t). This is a particular instance of what is referred to as duality between 
two stochastic processes in the probabilistic literature [13]. Thus in particular: 

and so one can have {ipnitY) < 0, while keeping varp„(i) = {{Pn{ty)) — {{Pnit)))^ — 
{Vn{ty) + {fn{t)) — {'Pn{'t)y ^0. It is also commouly accepted that the complex- 
valued nature of the trajectories of the field v^„(t) is needed in order to account for the 
importance of fiuctuation effects in low spatial dimensions, resulting in a slower decay 
for the total density of particles than what the naive law of mass action would predict 
(viz. r''/^ i/i in d < 2 for the reaction (1.7) [10, 14]). 

Henceforth, along the lines of [7, 8], we focus onto the single-site problem associated 
with (1.8), neglecting any spatial dependence. This simplification will allow us to better 
understand the role of the excursions of the field (p{t) in the complex plane. Setting 
K = 0, and absorbing the reaction rate A into the time scale, one ends up with the 
following Langevin-Ito equation for a single complex stochastic variable (p{t): 

^^-^(tr+Mtut), (1.11) 

where r]{t) is a normalized Gaussian white noise (see (1-3)). The initial condition (p{0) 
is real and non-negative (e.g., (p{0) = Pq if one starts from a Poisson distribution with 
density po for the original particle system). The conjugation symmetry (/?—>■ <^ ensures 
that the imaginary part of (p{t) averages over to zero, so as to maintain the reality 
and the non-negativity of {(p{t)) — {{p{t))). It has already been noticed by several 
authors [6, 7, 8] that (1.11) becomes linear in the variable Z{t) — l/(p{t). One thus 
obtains the explicit solution 

Z(t) = Z(0) e-*/2-'^(*) + /* e-{t-s)/2+ms)-B(t)] /i_;^2) 

Jo 

Rescaling time, and using the scaling property B{at) = y/aB{t) of Brownian motion, 
we obtain by identifying (1.5) and (1-12) the following identity between the stationary 
solution Z = hmt_>.oo Z{t) and the functional Z ior g — -\/2, i.e., D — 1 [8]: 

Z = lim Z(t) = 2Z\^. = 2 r e-'+'^^^'Mt. (1.13) 

t-5>oo ^ ^ Jo 

The above representation of the stationary solution has striking consequences. At 
the level of the original process (1.7), the stationary state is rather featureless, as there 
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just remains either zero or one particle, depending on the parity of the initial condition. 
By (1.10), and the corresponding equations for higher-order moments, this implies that 
= for any integer p > 2. Owing to (1.13), this reads {Z''^)\^^-^ = 0. The 
full distribution of Z will however turn out to be highly non-trivial. In particular, in 
stark contrast to what intuition backing up (1.8) or (1.11) could let us foresee, in the 
stationary state the variable (p — 2/Z never reaches the value (p — 0, characteristic of 
the absorbing state. Indeed, as already announced, Z lies within the unit disk. 

The setup of this article is the following. In Section 2 we present a self-contained 
investigation of the real exponential functional X (see (1-4)). The emphasis is put 
on the complementarity between the usual Langevin approach and another approach 
based on the similarity with the random recursions met in the study of Kesten variables 
and one-dimensional disordered systems. The main section of the paper (Section 3) 
is devoted to a detailed study of the imaginary exponential functional Z (sec (1.5)). 
We present numerical illustrations of the distribution of Z and investigate many facets 
of the problem by analytical means, including the relationship between the Langevin 
and Kesten approaches, the moments {Z^Z^), the weak-disorder and strong-disorder 
regimes, and the asymptotic behavior of the distribution near the unit circle. Section 4 
contains a brief summary of our findings. 

2. A warming up: the real functional 

This section is to a large extent intended as a warming up. It is devoted to a self- 
contained study of the real exponential functional (see (1.4)) 



The positive random variable thus defined is one of the exponential functionals of 
Brownian motion which have been investigated in probability theory, chiefly by Yor 
and his collaborators (see [1] for a review). It also appears in the physics literature, in 
the context of one- dimensional disordered systems [2, 3, 4]. The coupling constant g 
measuring the strength of noise, or, equivalently, the diffusion constant D (see (1-6)), is 
the sole parameter entering the definition of X. 

2.1. Langevin approach 

A first approach to study the random variable X consists in using Langevin equations. 
At this point it is useful to recall some elements of stochastic calculus [11, 15, 16]. 
A Langevin equation of the form 



is ambiguous as soon as it is non-linear, in the sense that the noise 77 (t) multiplies a 
non-trivial function h{X{t)) of the position X{t). This ambiguity due to the usage of 
a continuous-time formalism can be lifted in many ways. The two most useful and 
well-known prescriptions are the following (see [11, 16] for a detailed exposition): 




(2.1) 



dXjt) 



a{X{t)) + b{X{t)Ut) 



(2.2) 
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• Stratonovich prescription. The Langevin-Stratonovich differential equation 

[S] ^ = as{X{t)) + b{X{t)Ut) (2.3) 

can be essentially thought of as an ordinary differential equation. It is amenable 
to non-linear changes of variable according to the usual rules of differential and 
integral calculus. The main disadvantage is that X{t) and 77 (t) at the same time t 
are not independent. 

• ltd prescription. The Langcvin-lto differential equation 

[I] ^^a,{X{t)) + b{X{t)Ut) (2.4) 

has the advantage that the process X{t) and the noise r]{t) at the same time t are 
independent, so that one has e.g. {ai{X)) = in the stationary state. Equation (2.4) 
also provides a natural discretization of the process X{t). Considering discrete times 
t — ne so that Xn = X{tn), we obtain the recursion 

= X„ + ai{Xr,)e + 6(X„)C„+i, (2.5) 

where Cn+i is a Gaussian random variable, independent of Xn, such that (Cn+i) ~ 
and (Cn+i) — This discrete scheme can be efficiently used in a numerical 
simulation. The main disadvantage of the Ito prescription is that care must be 
exercised when making non-linear changes of variable. 

Both Langevin equations (2.3) and (2.4) describe the same stochastic process 
{X{t)} if their drift terms are related to each other by the correspondence formula 

a,iX) - as{X) = 1 A bixr = \ h{X)^. (2.6) 

The corresponding time-dependent probability density P{x,t) — {5{X{t) — x)) obeys 
the Fokker- Planck equation 

= -|M^)4|(''w|(k.)^))- (2.7) 

The stationary density of the process, fx{x) = limt^oo P{x, t), reads 

where xo is an arbitrary initial point, and A^i and Ns are normalization constants. 

It is now time to return to our functional X (see (1-4)). The Langevin-Stratonovich 
equation (1.2), i.e., 

[S] ^^l-X{t)+gX{tUt), (2.9) 
is equivalent to the Langevin-Ito equation 

[I] ^ ^1 + {D- l)X{t) + gX{t)r,{t). (2.10) 
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Indeed h{X) = gX and as{X) = 1-X yield ai{X) = 1 + {D - 1)X (see (1.6)). 

The Langevin-Stratonovich equation (2.9) can be readily integrated. We thus 
obtain the explicit stochastic representation of the process as 

X(t) = X(0) e-*+5^W + /* e-it-s)+9{B{t)-Bis)) (2.11) 

Jo 

In the i — > oo limit, the above expression loses the memory of its initial condition -'^(O). 
Moreover, using the stationarity of Brownian motion, the integral can be recast as an 
integral over t — t — s, which identifies with (2.1) or (1.4). We have thus shown that 
the exponential functional X represents the stationary solution of (2.9) or (2.10), i.e., 

X = \imt-,ooX{t). 

The Ito prescription allows us to directly read off the stationary mean value of X 
from (2.10): 

(X) = (2.12) 

This expression only makes sense for < 1. It is in agreement with the behavior of the 
time-dependent mean value {X{t)). Equation (2.10) indeed yields 
d{X{t)) 



dt 

whose solution, i.e.. 



^l + (D-l){X{t)), (2.13) 



{X{t)) = ^ + (X(0) - ^) e(^-^)*, (2.14) 

relaxes exponentially fast to (2.12) for D < 1, whereas it diverges for D > 1. The 
result (2.14) can be recovered by averaging (2.11) over the Brownian motion {B{s)}. 
Equation (2.8) leads to the following result for the distribution of the functional X: 

This expression [1, 3, 4] falls off exponentially fast as a; — )■ 0. It exhibits a fat tail 
at large x, as it decreases as a power law, with the continuously variable exponent 
— (1 + 1/Z^). As a consequence, the moment (X^) only converges for s < 1/D (see (2.30)). 
This provides another way to explain the divergence of (2.12) at D = 1. 

Finally, at par with the qualitative discussion given in the Introduction, one can 
interpret the emergence of a fat tail at large values of X and the possible divergence of 
the mean value {X) as consequences of the exchange of stability between the two fixed 
points of the dynamical system (2.9), i.e., the unstable deterministic one at X = 1 and 
a fluctuating one sitting formally at X = oo. 

Figure 1 shows a plot of the density fx for three values of D. 



2.2. Kesten approach 

An alternative approach to study the random variable X consists in using the similarity 
between the random integral (2.1) and random sums of products referred to as Kesten 
variables. Let us start with a reminder. A Kesten variable [17, 18, 19, 20, 21] is defined as 

^ = 1 + ^1 + 66 + ^166 + •••, (2.16) 
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Figure 1. Plot of the probability density fx of the real functional X (see (2.15)), for 
three values of the diffusion constant D. 



where the are i.i.d. positive random variables with probabihty density /^(C)- If the 
latter distribution is such that (In,^) < 0, the sum (2.16) is almost surely convergent, 
and it represents the stationary solution of the random recursion 

where is independent of Z„ and has density other words, we have the 

identity among random variables 

Z = 1 + ^Z', (2.18) 

where Z' is a copy of Z, and ^ is independent of Z'. The density fz{z) of the Kesten 
variable Z obeys the integral equation 

fz{z)= r^fdOfz('-^], (2.19) 



which cannot be solved in closed form in general. 

Let ,^min = a and ,^max = b he the smallest and largest values of ^ (lower and upper 
bounds of the support of f^), and similarly Zmm = A and Z^^x = B the smallest and 
largest values of Z. The condition (In^) < implies a < 1. We have then A = 1/(1 — a). 
If a < 6 < 1, we have B = 1/(1 — b) and fz has finite support. In the more interesting 
situation where a < 1 < 6, the support of fz extends to infinity. It is known that the 
distribution fz generically exhibits a fat tail, i.e., a power-law fall-off, of the form 

fz{z) ~ z-'^'+-\ (2.20) 

where the exponent a > is given by the condition [17, 19, 20] 

(D = 1- (2.21) 

The density fz{z) has been derived explicitly [18, 20] in the case where /^(O is a power 
law on an interval of the form [0, b] or [a, oo]. 
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The exponential functional X, defined in (2.1) as an integral of an exponential, 
can be viewed as a continuous analogue of the Kesten variable Z, defined in (2.16) 
as a sum of products. Similarly, the Langevin equation (2.9) is the continuous time 
analogue of the recursion (2.17). This correspondence, already referred to in [2, 3], can 
be understood quantitatively in the following way. Let us introduce a small time step e 
and split the integral in (2.1) as X = X^^'> + X^'^\ where 

= e-*+^^Wdt, = 2 e-*+^^Wdt. (2.22) 

To first order in e, we have X^^'> — e. Furthermore, setting t — e + s, we have 
B{t) = B{e) + B{s), so that X^^) = e'^+'^X', where X' is a copy of the variable X 
and ( = gB{e) is a Gaussian variable, independent of X', such that {() = and 
(C^) = 9^^ = 2-De. Putting everything together, the exponential functional X appears 
(in the s ^ limit) to obey the identity 

X = e + ^X', (2.23) 

with 

e = e-^+^, (2.24) 

where ( is Gaussian, such that (C) = and (C^) = 2De. Up to an unimportant global 
factor e, the exponential functional X therefore identifies (in the £ — > limit) with the 
Kesten variable Z generated by the random input variables ^ given by (2.24). This is 
precisely the Kesten variable investigated in [21], where the main emphasis is already 
put on the e ^ limit, and where the distribution (2.15) is derived in this limit. 

The above correspondence directly yields the fall-off exponent of the density of X. 
We have indeed (e*^ = e^'«')/2 = e^^'^ so that 

) = e*(^*-i)^ (2.25) 

The condition (2.21) thus predicts a = 1/D, in agreement with (2.15). 

The full distribution of X can actually be derived from the identity (2.23). Gonsider 
indeed the moment function 

M{s) = {X'). (2.26) 

Equations (2.23) and (2.25) can be respectively expanded in the £ — > hmit to yield 

M{s) = r)M(s) + s{C-')M{s -!)£ + ••• (2.27) 

and 

{C) = l + s{Ds~l)e + ---, (2.28) 

where the dots stand for terms of higher order in e. We are left (in the £ ^ limit) 
with the functional equation 

M{s-1) = {1- Ds)M{s), (2.29) 

whose normalized solution reads 

^(^) = T(^^y^ (Re.<l/i?). (2.30) 
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Setting s = 1, we recover the expression (2.12) for (X), provided D < 1. The density 
of X is given by the inverse MeUin transform 

fx{x)^ j^^^"~'M{s). (2.31) 

The expression (2.15) is recovered by summing the contributions of the poles of the 
integrand at s = 1/D + n for n = 0, 1, . . . The leftmost pole at s = 1/D is responsible 
for the power-law tail with exponent —{1 + 1/D). 



3. The imaginary functional 

We now turn to the main object of this paper, namely the distribution of the random 
variable Z defined by the integral (1.5), i.e., 

Z= / e-*+^^^(*)di. (3.1) 
Jo 

Prom a formal viewpoint, the complex variable Z can be viewed as the analytical 
continuation as g ^ \g of its real counterpart X, investigated in Section 2. This 
continuation amounts to changing the sign of the diffusion constant {D — > —D). It 
is however worth emphasizing that the study of Z is far more difficult than that of X, 
as most of the usual tools which are fit to investigate real random variables cease to 
work in the case of a complex random variable. 

Setting Z = X + iV, we are equivalently interested in the joint distribution f[x, y) 
of the two correlated real random variables 

e-^ cos{gB{t))dt, e'^ sm{gB{t))dt. (3.2) 

Jo Jo 

Let us start with a few general facts. The complex random variable Z lives inside the 

unit disk, since 

oo 

t 



\Z\ < / e-'dt = 1. (3.3) 
Jo 

Let us anticipate that its distribution has a smooth density f{x,y), whose support is 

the whole unit disk. Before we turn to more specific features, it is illustrative to first 
evaluate the mean (Z). Using the identity 

we obtain at once 

(Z) = r e-*(e^^^W)d^ = r e"^'''"'^* dt, (3.5) 
Jo Jo 



so that 



(Z) = (3.6) 
{X) = {Y} = 0. (3.7) 
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The above results obey the conjugation symmetry: Z = X+iY and Z = X — iY have the 
same law. The result (3.6) can be recovered by performing the analytical continuation 
D — )■ —D on (2.12). At variance with the latter result, the expression (3.6) depends 
smoothly on the diffusion constant D, decreasing from (Z) = 1 in the D ^ limit to 
{Z) ^ in the D ^ oo limit. These limiting regimes will be respectively investigated 
in Sections 3.5 and 3.6, whereas the main features of the dependence of the distribution 
oi Z on D will be studied in Section 3.7. 



3.1. Langevin approach 

It can be shown along the lines of Section 2.1 that the random variable Z = lim^^oo Z{t) 
represents the stationary solution of the Langevin-Stratonovich equation 

[S] ^^l-Z + igZrj{t), (3.8) 
where r]{t) is again a Gaussian white noise. The solution to (3.8) reads 

Z(t) = Z{0) e-*+i^?^W + /* e-(t-s)+i9{B{t)-B(s)) .3_9^ 

Jo 

A comparison with (1.12) shows that the identity (1.13) extends to finite times as 
Z{t) = 2Z(t/2) (again with D = 1). 

The Langevin-Ito equation corresponding to (3.8) is 

[I] ^^i-{D + l)Z + igZrj{t). (3.10) 

This equation yields (1 — (-D + l)Z) = in the stationary state. We thus readily recover 
the expression (3.6) of {Z). Equation (3.10) also provides an efficient discrete scheme. 
Considering discrete times t — ne, and adapting the recursion (2.5) to the present 
situation, we get 

Zn+l — S + Cn+lZn, (3-11) 

with 

in^l-{D + l)e + iCn, (3.12) 

and where Cn is a Gaussian random variable, such that (Cn) = and (C^) = 2Ds. 

Figure 2 shows a typical trajectory {Z{t)} starting from the origin, of length t = 25, 
for D = 1, i.e., the case of the binary annihilation process (1.7). The trajectory is 
generated by means of the scheme (3.11), with e = 10^'^. The main advantage of using 
the complex coordinate Z{t) (instead of the original field variable ip{t) [7, 8]) is that 
unbounded excursions are avoided, as we have \Z{t)\ < 1. The rest of this section is 
devoted to the distribution of the functional Z, i.e., equivalently, of the generic point of 
such a trajectory in the long-time regime. 
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Figure 2. A typical trajectory {Z{ty\ of length t ^ 25 in the unit disk, for D = 1, 
i.e., the case of the binary annihilation process (1.7). 



3.2. Kesten approach 

The random variable Z can be alternatively described as a complex Kesten variable. 
Along the lines of Section 2.2, it can indeed be shown that Z obeys the identity (in the 
e limit) 

Z = e + ^Z\ (3.13) 

where Z' is a copy of Z, whereas 

e = e-^+^S (3.14) 

and C is Gaussian, such that {Q = and {(^) = 2De. In other words, Z represents (in 
the e — )■ limit) the stationary solution of the random recursion 

(3.15) 

with 

^„ = e-^+'f". (3.16) 

The variable Z therefore identifies (up to an unimportant global factor £, and in the 
£ — )■ limit) with the Kesten variable generated by the complex input variables ^n- 

It is worth noticing the close resemblance between the random recursions (3.11) 
and (3.15), respectively corresponding to the Langevin and Kesten approaches. Both 
discrete schemes become equivalent to leading order in the e — )■ limit. The Langevin 
scheme involving ^„ appears in this regime as a suitably linearized form of the non- 
linear Kesten scheme involving The Langevin scheme is more suitable for extensive 
numerical simulations, as it does not involve the complex exponential function. 
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3.3. Moments 

In the case of the real exponential functional X, the explicit form (2.15) of the 
distribution fx could be derived either from the Langevin approach, as the stationary 
solution (2.8) of the Fokker-Planck equation, or from the Kesten approach, by means of 
the identity (2.23) in the e — )■ limit. In the present case of the imaginary exponential 
functional Z, neither route leads to an explicit expression for the density f{x,y). 

The Kesten approach however directly yields an interesting piece of information, in 
the form of a recursion relation for the moments (Z^Z^). We recall that the distribution 
of a real random variable X is entirely characterized by the single array of moments (X^) 
{k — 0,1, . . .) (leaving aside questions related to convergence), whereas a double array 
of moments I^Z^Z^) (or, equivalently, {X^Y^y) is needed to characterize the distribution 
of a complex random variable Z — X ^ SY . 

The analysis of the moments follows the line of thought of Section 2.2. For 
A; = 0, 1, . . . integer, we expand the identity (3.13) as 

Z^ = i^{Z'f + ki^-^{Z'f-^e + • • • (3.17) 

We thus obtain an expansion similar to (2.27), i.e., 

{z'^z') = {e-^){^z'z') 

+ (k{i''-^l'){Z''-^Z') + + • • • (3.18) 

Furthermore, we have 

^l-{k + l + {k-lfD)e + --- (3.19) 

Inserting the latter expansion into (3.18), we are left (in the £ — > hmit) with the 
recursion 

{k + l + {k- ifD) {Z^Z^) = k{Z^-^Z^) + 1{Z^Z^-^) . (3.20) 

The formula (3.20) is one of the key results of this work. It allows one to evaluate 
the moments (Z'^Z') in a recursive way. All of them are rational functions of the diffusion 
constant D. They are regular over the whole range of positive values of D. Owing to 
the symmetry between Z and Z, we have [Z^Z^) — [Z^Z^). so that we can restrict 
ourselves to A; > L The expressions of the first moments (for k + I < &) are listed in 
Table 1. 

The moments {Z^) can be evaluated in closed form. Equation (3.20) indeed boils 
down for / = to the recursion 

{kD + l){Z^) = {Z^-^), (3.21) 

which can be viewed as the analytical continuation of the functional equation (2.29). 
We thus obtain 

^ '~ I, ~ D^T{k + 1 + 1/ D)- ^"^-^^^ 
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moment 


expression 




\ 

D+l 




1 

(D+1)(2Z)+1) 


{ZZ) 


1 

D+l 


1 '73\ 

kz^z) 


1 

(D+1){2D+I)(3r>+1) 

4Z)+3 
(D+1)(D+3)(2D+1) 


\Z-'z) 

{Z^Z^) 


1 

(D+1)(2D+1)(3D+1)(4D+1) 

(£)+l)2(D+3){2D+l)(3D+l) 
4_D+3 
(D+l)(D+3)(2£l+l) 


{Z^Z) 
\Z z ) 


1 

(D+l)(2D+lH3ri+l)flD+l)(5r)+l) 

lHiJ- + 13«+3 
(D+l)2(D+3)(2Ll+l)(3D+l)(4£l+l) 

36D3 +93D2 +68D+ 15 
(D+l)2(D+3)(D+5)(2£)+l)(3D+l) 


{Z^'Z) 
{Z^Z^) 
{Z^Z^) 


1 

(D+1)(2D+1)(3D+1)(4D+1)(5D+1)(6I3+1) 

25Z)2 + 16_D+3 
(D+l)2(D+3)(2i?+l)(3D+l)(4D+l)(5D+l) 
288D3+544D2+279D+45 


(D+1)(D+3)(D+5)(2D+1)(2D+3)(3D+1)(4D+1) 
361)3 +931)2 _,_gg£)_,_;^5 

(D+l)2(D+3)(D+5)(2£)+l)(3D+l) 



Table 1. Expressions of the moments {Z'^Z^) iovk+l < 6 as a function of the diffusion 
constant D. 



For D = 1, the above expression simplifies to 

This result also holds for negative k. Setting k = — p, with p > 2 integer, we obtain 
{Z''P)\jj^^ = 0. We thus recover a result mentioned in the Introduction, concerning the 
null factorial moments for the number of particles in the steady state for the binary 
annihilation process (1.7). 

The joint moments [X^Y^) of the real variables X and Y can be obtained as linear 
combinations of the moments {Z'^Z''), using X = {Z + Z)/2 and Y = i{Z — Z)/2. The 
conjugation symmetry ensures that f{x,y) = f{x, —y), so that {X^Y^) = whenever / 
is odd. The expressions of the first non-vanishing moments {X'^Y'') (for / even and 
A; + / < 6) are listed in Table 2. 



3.4- Fokker- Planck equation 

The next natural step in our analysis consists in writing down the Fokker-Planck 
equation obeyed by the distribution of the variable Z, or, equivalently, by the joint 
distribution of the real variables X and Y. We shall successively do so using various 
representations. 
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moment 


expression 


/ '\r\ 


I 

D+l 




1 

2D+1 

D 

(D+1)(2D+1) 


{X ) 


9r)^+ior>+3 

(D+1)(D+3)(2D+1)(3D+1) 
3D 

(D+3)(2D+1)(3Z)+1) 


{X Y ) 


18Z)*+60D'^+58D2+22r>+3 
(D+1)2(D+3)(2D+1)(3Z)+1)(4D+1) 

D(6Z)2+8D+3) 
(D+1)(D+3)(2D+1)(3D+1)(4D+1) 

3Z)2(6D2+8D+3) 
(D+l)2(D+3)(2D+l)(3i3+l)(4D+l) 


{XY } 


450i:i*+940D3+606D2+158D+15 

(D+1HD+3)(D+5)(2D+1)(3D+1)(1D+1)(5D+1) 

(D+l)2(£l+3)(£>+5)(2L>+l)(3I3+l)(4D+l)(5D+l) 
3D2(3od3+86D2+63D+15) 

(D+1)2(D+3)(D+5)(2D+1)(3D+1)(4I?+1)(5I?+1) 


(X^Y^) 


2700Z)7+16740ri*'+35940ri-"''+36182D4 + 19056D-H5383Z)2+774£)+45 


(D+1)2(D+3)(D+5)(2D+1)(2D+3)(3D+1)(4D+1)(5Z)+1)(6D+1) 

3D(180D5+756D'*+1044Z)3+584Z)2_|_i48D+15) 
(D+1)(D+3)(D+5)(2D+1)(2D+3)(3Z)+1)(4Z)+1)(5Z)+1)(6D+1) 

3D2(180D'*+576D^+638Z)2+279D+45) 
(£l+l)(Z)+3)(D+5)(2r>+l)(2D+3)(3D+l)(4D+l)(5D+l)(6D+l) 

15D3 (180D''+576D3 +6381)2 +279D+45) 
(D+1)2(D+3)(D+5)(2D+1)(2£>+3)(3D+1)(4I?+1)(5D+1)(6D+1) 



Table 2. Expressions of the moments [X'^Y'') for I even and A; + Z < 6 as a function 
of the diffusion constant D. 



We start by considering the two-variable characteristic function 

q) = (e^^+^^) = £ {Z'ZY^. (3.24) 

The recursion (3.20) is equivalent to the partial differential equation 

/ d dV d d \ 

with $(0,0) = 1. For q — the above equation boils down to the solvable ordinary 
differential equation 

Dp(t>"{p) + {D + 1)0' (p) -0 = (3.26) 

for the characteristic function (j){p) = (e^^), where accents denote differentiations. We 
thus obtain 

<P{p) = r(l + l/D) (I) ''^^'""^ lyn (2^1) , (3.27) 

where Ii/d is the modified Bessel function. This result agrees with the expression (3.22) 
of the moments {Z''). 
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The characteristic function "^{u^v) = (^e'^^'^'"^) is equal to g), with p — 
{u — iv)/2 and q— {u-\-iv)/2. It therefore obeys the partial differential equation 

" ("I - "I;)' - "I; - + ") *("• = °- 

with ^(0, 0) = 1. The latter equation is equivalent to the following equation for the 
joint density of X,Y: 

The above equation can be recovered by noticing that /(x, y) is the stationary solution of 
the two-dimensional Fokker-Planck equation associated with the two coupled Langevin- 
Stratonovich equations 

dX{t) 



dt 
dY{t) 



^l-X{t)-gY{t)rj{t), 

-Y{t)+gX{t)r]{t). (3.30) 



dt 

A single real noise ri{t) enters both equations (3.30), and so the above system of Langevin 
equations is degenerate. As a consequence, the second-order differential operators which 
appear in either (3.25), (3.28), or (3.29) are of parabohc nature. This phenomenon is best 
exhibited in yet another form, using the polar coordinates Z — R&^. The system (3.30) 
then becomes 

dR{t) _ 



dt 

de{t) sme{t) 



R{t) + COS e{t), 

+ 9V{t): (3.31) 



dt R{t) 

showing explicitly that the noise only acts in the angular direction. This property was 
already noticed in [7, 8] in the context of binary annihilation (see (1.7)), i.e., for D — 1. 
Coming back to the stationary state, still denoting Z — R&^, and setting 

P=fe-^ g=|e^ (3.32) 

the characteristic function ^{p, q) becomes the generating function 

n{li,u) ^ {e"^-"), (3.33) 

where 

X^ = ^{Ze-''' + Ze''') =Rcos{e-v) = Xcosi/ + F sini/ (3.34) 

is the projection of Z onto an axis at an angle u. The function f2(/x, u) obeys 
( d \ 

I D— - fi— + iicosuj ^0, (3.35) 

with f2(0, z/) = 1 for all i^. 

The partial differential equation (3.35) is another central result of this work. It 
will be used in the following in order to investigate the density of Z in various limits. 
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We shall discuss the weak- disorder regime {D — >■ 0) in section 3.5, the strong-disorder 
limit {D — >■ oo) in sections 3.6 and 3.9, and the asymptotic behavior of the density of Z 
near the imit circle in section 3.8. In spite of its relatively simple form, we have not 
been able to solve the parabolic partial differential equation (3.35) in full generality. We 
nevertheless find it worth to sketch below two approaches we have tried. 

Fourier-series expansion. The generating function Q(/i, u) is 27r-periodic in the angular 
variable u. It is henceforth natural to expand it as a Fourier series: 

oo 

n{fx,u)= ^n(/^)e'"^ (3.36) 

n=— oo 

The Fourier coefficients have the symmetry u!_n{n) = cOnif^)- They obey the coupled 
differential equations 

^Dn^ + 1^-^ UJnil^) = I {l^n+l{l^) + UJn-l{l^)) , (3.37) 

with a;„(0) = 5n,o- For D — 0, (3.37) coincides with the differential recursion relation 
obeyed by the modified Bessel functions Ini/j)- Looking for asymptotics similar to those 
of Bessel functions, we obtain the following behavior of the amplitudes a;„(/x) at small 
and large values of for all n > 0: 

, , r(i + i/D) f Y , 

"-^^^" n!r(l+l + l/D) [2d) 

uMr^e^" ()">!)■ (3.38) 

A more advanced analysis of the large-// regime will be performed in Section 3.8. 

Expansion in trigonometric polynomials. The generating function 0(/x, v) can alterna- 
tively be expanded as a power series in of the form 

oo k 

n(l^,u)^J2ljPk(cosu). (3.39) 

A;=0 

Equation (3.33) yields 

Pfc(cosi/) = (X^) = {{Xcosu + Ysinu)'') = \{{Ze-''' + Ze''')''). (3.40) 

The Pk are polynomials of degree k in x = cos p. They satisfy the differential recursion 
relation 

D{{x^ - 1)P^' + xP^) + kPk = kxPk-i, (3.41) 
with Po{x) — 1. Alternatively, as a consequence of (3.40), we have 

P,(cosi.) = ^ E (^) {Z'-'Z^) cos((2j - k)u). (3.42) 
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Both (3.41) and (3.42) together with (3.20) or Table 1 consistently yield 

. , , , X , > X I D 

P,{X) = — — , P2ix) 



D + V ' (L' + 1)(2D + 1)' 

P (x) = (D + ^)^' + mD + l)^ .3 43. 

{D + l){D + 3){2D + l){3D+iy ^ ' 

and so on. 

The polynomials Pk{x) do not seem to belong to any of the known standard families 
of polynomials. In particular, they do not form a set of orthogonal polynomials. This 
can easily be checked using the fact that orthogonal polynomials must obey a three-term 
linear recursion [22]. 

An asymptotic analysis of the large-order behavior of the polynomials Pk{x) will 
also be done in Section 3.8. 

3.5. Weak-disorder regime (D ^ 0) 

In this regime, the complex exponential involved in (1.5) or (3.1) or, equivalently, the 
trigonometric functions involved in (3.2), can be expanded in powers oi g = y/2D. We 
thus obtain 

Z^l + i^/DU-DV + ---, (3.44) 

i.e., 

X = 1-DV + ---, Y = ^/DU + ■■■, (3.45) 
where U and V are the following random integrals: 

poo roo 

U^V2 e-*B(t)dt, e-*B(t?dt. (3.46) 

JO Jo 

The variable U is Gaussian, and therefore entirely characterized by its second moment 
(t/^) = 1. This is however not the end of the story. The variables U and V indeed have 
a non-trivial joint distribution supported by the parabolic domain V > As these 

variables are respectively defined as a linear and a quadratic functional of Brownian 
motion, one way to investigate their joint distribution consists in evaluating the 
corresponding characteristic function as a Gaussian functional integral. We prefer to use 
a more direct approach, based on the solution of the partial differential equation (3.35) 
in the relevant regime, i.e., D 1 and u <^1. We arc thus led to consider a simplified 
version of (3.35), where the cosine function has been expanded to second order: 

The latter equation can be solved for arbitrary D if, inspired by the property that U is 
Gaussian, we make a Gaussian Ansatz of the form 

^a{n)e-^^>'^''\ (3.48) 
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We are left with two ordinary differential equations for a(//) and namely 

with a(0) = 1 and 6(0) = 0. The equation for 6(/x) is a Riccati equation, and is therefore 
solvable. Indeed, setting 

we obtain for ■0(a*) the second-order linear differential equation 

'2D 1 
H An 



r{i^) = [ — -—2]^(iA. (3.51) 



whose regular solution is 

V^(^) = {2^2Di^ . (3.52) 

Skipping details, we obtain after some algebra 



^Jlo{2^/2Djl) Io{2^2Dn) 

where Iq and Ji are modified Bessel functions. In the following we use various properties 
of Bessel functions (differential equations, series and asymptotic expansions, location of 
the zeros) exposed e.g. in [23]. 

Now, inserting (3.45) into (3.33) and (3.34), and expanding the result for D and u 
small, we obtain 

fl{ii, v) = ^exp + VD/iuU - D/iV - ^iiu^ H ^ ^ (3.54) 



In the limit where /i is large, at fixed values of the combinations p — —yDfj^u and 
(7 = D/i, the leading term fj, of the above expansion is divergent, whereas all the other 
terms remain finite, and the higher-order terms which are not written down vanish. As 
a consequence, a direct identification with (3.48), (3.53) yields an explicit expression for 
the characteristic function 



F(p, a) = {e-^^-^^) = , ^ exp [ ^''^^^^ _ (3 55^ 

4U2^) '^\2aIo{2V2^)' J 

This expression is another key result of this work. It contains the full joint distribution 
of the rescaled variables U and V. Table 3 gives the rational values of the first few 
moments {U^V^) (for k even and k,l < 4). The positive variables C/^ and V appear as 
strongly positively correlated. 

The Gaussian marginal distribution of U is recovered as F(p, 0) = {e~f^) = e^^^^. 
The marginal distribution of the quadratic functional V is given by the inverse Laplace 
transform of F(0,cr), i.e., 

^ . . r da e"" 
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moment 


value 


moment 


value 


moment 


value 


{V) 


1 
2 

20 

3 

95 
3 


(f/V) 


1 

7 

3 
25 

3 
613 
15 
2305 

9 




3 
11 

164 
3 

1726 
5 

39751 
15 



Table 3. Values of the moments {U^V^) for k even and k,l < 4. 



The behavior of this density at small v is given by a saddle point at a large positive 
value (T ~ 1/ (2v'^) of the integration variable. We thus obtain the exponential fall-off 

~e-i/(2.) (^-^0). (3.57) 

Its behavior at large v is determined by the rightmost singularity of the integrand 
in (3.56), which is a branch-cut singularity at a = —<Ji, with o"i = jf/8 ~ 0.722897, 
where ji ~ 2.404825 is the first zero of the Bessel function Jq. We thus obtain the 
exponential fall-off 

/y(i;) -e-"^'' (v^oo). (3.58) 

Figure 3 shows a plot of the density /y(f ), obtained by means of a numerical integration 
of the formula (3.56) along the imaginary a axis. 

1.5 




V 

Figure 3. Plot of the probability density fv of the scaling variable V. 
3.6. Strong- disorder regime (D ^ oo) 

In this regime, the integrand in the expression (1.5) or (3.1) involves a large and hence 
rapidly varying phase. The variable Z can therefore be expected to behave as a complex 
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Gaussian variable with zero mean and a small width. Still denoting Z — Re^^, we have 
indeed {R^) = {ZZ) (see Table 1). 

The convergence of Z to a complex Gaussian variable can be studied by means 
of the moments investigated in Section 3.3. The form of the recursion (3.20) implies 
that the moments behave as {Z^Z^) ~ x^£)max(fc,o large D. As a consequence, the 
diagonal moments {R^^) = {{ZZY) are the leading ones at large D, whereas the 

non-diagonal ones {Z'^{ZZY) ~ [n — 1,2,.. .) are more and more subleading. 

More precisely, setting 

{Z-{ZZf) = (3.59) 
the recursion (3.20) is equivalent to 

{Dn^ + 2k + n)ai") = D{k + n)af-^^ + fcaK'^ {n > 1). (3.60) 
To leading order as D — )■ cxd, the latter equation boils down to 

n'at^^{k + n)at'\ (3.61) 
so that the amplitudes have the asymptotic values 

We thus obtain in particular 

{{ZZ)'') ^ (3.63) 

We read off from this result that the positive variable S = R^ = ZZ = + Y'^ is 
asymptotically exponentially distributed at large D, with density fs{s) ~ De^^'^. In 
other words, Z is asymptotically a complex Gaussian variable, with the isotropic density 

/(x,y)«-e-^«^ (3.64) 

TT 

The analysis of the recursion (3.60) can be pursued to derive a systematic 1/D- 
expansion beyond the asymptotic values (3.62). For the diagonal moments, we thus 
obtain after some algebra 

(3,5) 

3. 7. Qualitative features 

In this section we make a pause in our analytical investigations and turn to a qualitative 
discussion of various features of the distribution of the variable Z, supported by 
numerical results. 

The global picture is provided by Figure 4, showing color level plots of the density 
f{x,y) of the variable Z = X + iY in the unit disk, for several values of D. The 
data presented here and throughout the following have been obtained by numerically 
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generating very long time series of the discrete process (3.11), with typically e = 10~^ 
and n = 10^°. First of all, as the distribution is symmetric w.r.t. the X-axis, its unique 
maximum is on the real axis. Figure 5 shows a plot of the observed location Z^ax of this 
maximum (black curve) and of the mean (Z) (red line, see (3.6)), against D/{D + 1). 
Both quantities qualitatively follow the same pattern, decreasing from 1 in the D — ?■ 
limit to in the D ^ oo limit. 



D=l/4 D=l/2 D=l 




Figure 4. Level plots of the density f{x,y) of the variable Z = X + iY in the unit 
disk, for several values of D. 

The full shape of the distribution is observed to exhibit a smooth dependence on 
the diffusion constant D, continuously interpolating between the weak-disorder and 
strong-disorder regimes, respectively studied in Sections 3.5 and 3.6. 

• For small D, the distribution of Z is concentrated in a very elongated domain 
which is stretched along the unit circle and fits its curvature. This very anisotropic 
distribution is best illustrated on the first panel of Figure 4, corresponding to 
D = 1/4. The picture corroborates the fact that Y has a variance {Y"^) ^ D, 
and therefore a vertical extent scaling as \/D, whereas the horizontal extent of 
1 — X ^ DV is much smaller. The curvature reflects the property that f/^ and V 
are strongly positively correlated. 

• For large D, Z is predicted to be asymptotically a complex Gaussian variable such 
that (-R^) ^ 1/D. Its distribution is therefore essentially concentrated in a small 
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V 
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N 



0.2 



<z> 




0.4 0.6 
D/(D+1) 

Figure 5. Plot of the location Z^aax of the maximum of the distribution of Z and of 
the mean {Z) (see (3.6)), against D/{D + 1). 



and nearly circular domain near the origin. This nearly isotropic distribution is 
illustrated on the last panel of Figure 4, corresponding to D = 8. 

The expressions of the variances of X and Y , 

""^^=(D + 1)^(2D + 1)' ""^^=(^ + 1)^2/5 + 1)' ^'-''^ 

are readily derived from Table 1. Both variances go to zero at both endpoints (-D — )■ 
and D — )■ oo). For small varX D"^ vanishes much faster than varF ^ D, 
reflecting the very anisotropic nature of the distribution. To the contrary, at large D 
both variances are asymptotically equal to each other, as varX varF ^ 1/{2D), 
reflecting the nearly isotropic nature of the distribution. The variances are maximal for 
intermediate values of the diffusion constant. The maximum (varX)max = (5-\/5 — 
ll)/2 ^ 0.090169 of varX is reached for D = {I + v^)/2 ^ 1.618033, whereas 
(varF)^ax = 3 - 2^2 ^ 0.171572 is reached for D = 1/^2 ^ 0.707106. It is worth 
noticing that varF is invariant if D is changed to its dual D, such that 

DD = \. (3.67) 

The maximum (varF)inax is consistently reached at the fixed point D = l/\/2 of this 
duality transformation. Figure 6 shows a plot of both variances against D/(D + 1). 

The features discussed so far are corroborated by the shape of the marginal 
distributions fx and fy of the real variables X and Y. These distributions, obtained 
as numerical histograms, are shown in Figure 7 for several values of D, chosen so as 
to form dual pairs (see (3.67)). The dependence of the density of X on the diffusion 
constant is observed to evolve smoothly from a narrow distribution around X = 1 at 
small D, whose limiting shape is given by the density fv of the rescaled variable V, to 
a narrow Gaussian near X = at large D, whose variance scales as varX ^ 1/{2D). 
The distribution of Y behaves differently. It indeed becomes a narrow Gaussian both at 
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Figure 6. Plot of the variances varX and varl^ (see (3.66)) against D/{D + 1). 



weak and at strong disorder, whose variance respectively falls off as var Y D {D ^ Q) 
and vary ^ 1/{2D) [D — t- oo). The densities of Y for pairs of dual values D and D of 
the diffusion constant are observed to be strikingly close to each other. 




Figure 7. Plot of the densities fx (left) and fy (right) of the marginal distributions 
of X and F, for values of D chosen so as to form the dual pairs (8, 1/16), (4, 1/8), (2, 
1/4), and (1, 1/2) (see (3.67)). 



3.8. Asymptotic behavior near the boundary: an essential singularity 

We now resume our analytical investigations and turn to more advanced features, namely 
the asymptotic analysis of the distribution of Z in two regimes where it is exponentially 
small: near the boundary (unit circle) for all values of D (in this section), and for 
large D at any point inside the unit disk (in Section 3.9). 

For the time being, our goal is to investigate the distribution of Z near the unit 
circle. More precisely, we are interested in the behavior of the density fij{Xy) of the 
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variable defined in (3.34) as the unit circle is approached from inside {X^ — >■ 1), for 
a fixed value of the angle v. The definition (3.33) implies 

f.{X.)= I ^e-'^^^n(/x,z/). (3.68) 

Wc are thus led to study the asymptotic behavior of the characteristic function fl{iJ,, v) 
at large and fixed v. 

It is worthwhile to first have a look at the simplified situation considered in 
Section 3.5, i.e., v <^ 1. In this regime, inserting (3.53) into (3.48), and expanding 
the result for large we obtain 



n(//, v) = exp |/x - ,j2Dii + ^ 



+ iln(327r2D/.)-^ + ---|. (3.69) 

Inspired by the above expansion, we look for a solution of the full partial differential 
equation (3.35) which behaves at large as 

Q(//,i/) =exp(//->l(z/)^+---). (3.70) 

Inserting this expansion into (3.35), we are left with the ordinary differential equation 

for the amplitude A(iy). Let us assume that the value A(0) = \J2D at z/ = 0, which can 
be read off from (3.69), holds for all values of D. This hypothesis will be tested against 
numerical results in a while. We thus obtain 

A(i/) = ^(^D + 2-2cos0 . (3.72) 

The fall-off of the distribution of Z near the unit circle can be estimated by 
inserting (3.70) into (3.68): 

For Xi, — >■ 1, this integral is dominated by a saddle point at // A(i/)^/(4(l — Xi,y). 
This yields the estimate 

i.e., 

MX.) - exp (- ^ 2D(l-xJ ) ■ ^'-''^ 

This result confirms that the distribution of Z is supported by the whole unit disk, 
and predicts that it vanishes exponentially fast near any point of the unit circle. The 
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behavior of the densities of the variables X and Y near their endpoints is given by the 
following simpler formulas, respectively corresponding to = 0, = ±7r, and v — ±7r/2: 

~exp(-— ^— ) {x^l), 




The first of these estimates matches the behavior (3.57) of the density of V. 

The dependence of the amplitude A{u) on D is not monotonic. It diverges 
at both endpoints (D — )■ and D oo), and therefore reaches a minimum, 
^min('^) = 4^1 — cos(z//2), for an angle-dependent value of the diffusion constant, 
-Dmin(^) = 2 — 2cos(z//2). This non-monotonic behavior holds for any angle except 
along the real axis {u = 0), where wc have ^(0) = \/2D. It is reminiscent of the 
behavior of varF discussed in Section 3.7. 

The essential singularity (3.75) of the distribution of Z near the unit circle is 
unobservable in numerical data. Nevertheless, it has an observable consequence. It 
indeed reflects itself in the large-order behavior of the polynomials Pk{cosu) — (X^) 
(see (3.40)). We have 

Pfc(cos u) = Xl f,[X,) dX,. (3.77) 

For large /c, the integral is dominated by either of the endpoints. For an angle v < 7r/2, 
the leading endpoint is X^, ^ 1, so that the integral can be estimated as 

P,(cos u) ^ £ X^, exp (-^0)^^ dX^. (3.78) 

The integral is dominated by a saddle point a,t X^ 1 — A{u) / {2y/k) , yielding 

Pfe(cos v) ~ exp (^-A{u)Vk) . (3.79) 

We therefore predict a stretched exponential decay of the polynomials Pfe(cos u), at any 
fixed angle u. Such a behavior can easily be observed numerically. The recursion (3.20) 
and the expansion (3.42) indeed allow a very accurate numerical evaluation of at least 
the first 1,000 polynomials. Figure 8 shows a plot of the amphtude A{i') against the 
reduced angle u/tt < 1/2, for several values of the diffusion constant D. The striking 
agreement between the prediction (3.72) for the amplitude (full lines) and the result of 
a quadratic fit of all the first 1,000 polynomials as \nPk(cosiy) = -Ay/k + Blnk + C 
(symbols) shows in particular that the hypothesis made above (3.72) is valid. 



3.9. Large deviations at large D 

We close our investigation by returning to the regime where the diffusion constant D is 
large. In this situation, it has been shown in Section 3.6 that the bulk of the distribution 
of Z is a narrow isotropic Gaussian, whose density is proportional to e^^^^ (see (3.64)). 




Figure 8. Plot of the amplitude A{i') against the reduced angle i'/tt, for three values 
of the diffusion constant D. Full lines: analytical prediction (3.72). Symbols: result of 
a fit of the first f ,000 polynomials Pfc(cosz^) (see text). 



Furthermore, the estimate (3.75) of the distribution near the unit circle also becomes 
isotropic at large D, and it falls off at e~^^^'^^^^^-^\ These observations suggest that the 
distribution follows an exponential large-deviation estimate of the form 

/(a;,i/)^/(i?)~e-^^(«) (3.80) 

for large D, all over the unit disk. 

Our goal is to show that such an estimate indeed holds, and to derive the large- 
deviation function S{R). First, if the distribution of Z is nearly isotropic, we have 

Q{fi,u) ^27r Rf{R)Io{fiR)dR, (3.81) 
Jo 

for all values of u. Then, if the distribution obeys the formula (3.80), the above integral 
can be evaluated for large fi (i.e., /i of the order of D) by means of the saddle-point 
approximation. Setting fi = Dx, we obtain 

fi(/i, z/) ~ e^^("\ (3.82) 

where the functions S{R) and S(x) are the Legendre transforms of each other: 

J:{x) + S{R) = xR, x = S'{R), R = E'{x). (3.83) 

The formulas (3.80) and (3.82) are just exponential estimates. In order to proceed, 
we have to consider the angular dependence of the pre-exponential factor in the latter 
estimate. Setting 

fi(/i,z/) ^i7(z/)e^^("), (3.84) 

the partial differential equation (3.35) yields the ordinary differential equation 

H"{u) + {x cos u - E)H{u) = (3.85) 

for the amplitude H{h') at fixed x = fj,/D, with E = xR = jiR/D. Equation (3.85) is 
known as the Mathieu equation [23]. The requirement that H{v) be periodic and positive 
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implies that E{x) is the ground-state eigenvalue of the latter equation. Introducing the 
Fourier series 

oo 

H{y)^ E ^ne''^^ (3.86) 

the Mathieu equation (3.85) is equivalent to the three-term recursion 

2{n^ + E)K^x{K+i + K-i), (3.87) 

which is the appropriately rescaled form of the recursion (3.37). The ground-state 
eigenvalue E{x) is characterized by the property that the above recursion has a positive 
normalizable solution. The function E{x) being known, S(x) is obtained by integrating 

the relation T!{x) = E{x)/x with E(0) = 0, and finally we have R = E{x)/x and 
S{R) = E{x) — S(x). We have thus obtained the large-deviation function S{R) in 
parametric form. Analytical expansions can be performed for x ^ and x ^ oo, 
respectively corresponding to i? — > and i? — )■ 1. 

For a; ^ 0, the amplitudes /i„ ~ a;" of the harmonics fall off very fast. Taking 
successive harmonics into account, we can derive the expansions 

^^''^ + + 

7 SQ^ 

5(i?) = i?^ + ^i?^+||i?^ + --- (i?^0). (3.88) 

For X — )■ oo, the function H{v) becomes peaked around z/ = 0. In this regime, (3.85) 
maps onto a weakly anharmonic quantum-mechanical oscillator. Treating the leading 
non-linearity in to first order, we obtain 

£;(a;) = X - + ^ H {x ^ oo), 

^^^^ = W^) + \ ^""^^ -R) + C + --- {R^l). (3.89) 

The leading behavior of S{R) at both endpoints matches the known results recalled 
in the beginning of this section. The coefficient 1/8 of the logarithmic term in S{R) 
near R = 1 matches the similar coefficient in (3.69). 

The above scheme is easily implemented numerically. This yields e.g. the constant 
C ~ —0.603. Figure 9 shows a plot of the function S{R) thus obtained (left), and of 
the small difference 5S{R) = S{R) - i?V(l - R^) (right). 



4. Discussion 



This work has been devoted to a detailed investigation of a random complex integral Z, 
which provides the most natural example of an imaginary exponential functional of 
Brownian motion. Both in the warming-up part (Section 2) devoted to the real 
functional X and in the main part (Section 3) devoted to Z, the main emphasis 
has been put on the complementarity between the more traditional approach, based 
on Langevin equations and stochastic calculus, and a more original one, based on 



On an imaginary exponential functional of Brownian motion 29 




R R 

Figure 9. Plot of the large-deviation function S{R) characterizing the distribution 

of Z at large D (left), and of the difference 5S{R) = S{R) - R^/{1 - R^) (right). 



discrete random recursions and Kesten variables. Even though neither of these routes 
leads to the solution of the problem, i.e., to an explicit expression for the probability 
distribution function of Z, we have gathered many partial results concerning various 
facets of this distribution in various regimes. Let us emphasize in particular the recursion 
relation (3.20) for the moments {Z'^Z'-), and the occurrence of an essential singularity 
of the form (3.75) of the asymptotic behavior of the density near the unit circle. 

Apart from its intrinsic interest, the imaginary exponential functional Z has close 
connections with the realm of reaction-diffusion processes. This fact, already underlined 
in [8], was our original motivation for pursuing the present work. The ubiquity of 
imaginary noise Langevin equations illustrates the potential pitfalls one faces when 
building up phenomenological stochastic equations in out-of-equilibrium situations. In 
fact, within the Doi-Peliti formalism for reaction-diffusion processes, it can be shown [24] 
that degenerate Langevin equations for a complex-valued field are the rule rather than 
the exception. The present study will hopefully trigger further interest in such stochastic 
processes. 
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